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ABSTRACT 

We investigate the equilibrium properties of self-gravitating magnetized clouds 
with polytropic equations of state with negative index n. In particular, we con- 
sider scale- free isopedic configurations that have constant dimensionless spherical 
mass-to-flux ratio and that may constitute "pivotal" states for subsequent dy- 
namical collapse to form groups or clusters of stars. For given T = I + 1/n, 
equilibria with smaller values of are more flattened, ranging from spherical 
configurations with A^ = oo to completely fiattened states for A^ = 1. For a 
given amount of support provided by the magnetic field as measured by the di- 
mensionless parameter Ho, equilibria with smaller values of F are more fiattened. 
However, logatropic (defined by F = 0) disks do not exist. The only possible 
scale-free isopedic equilibria with logatropic equation of state are spherical uni- 
formly magnetized clouds. 

Subject headings: ISM: clouds — stars: formation — MHD 
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1. INTRODUCTION 

The stage leading up to dynamic collapse of a magnetically subcritical cloud core to a 
protostar or a group of protostars is believed to be largely quasi-static, if the responsible 
process is ambipolar diffusion (e.g., Mestel & Spitzer 1956, [Nakano 1979| , [Lizano &: Shu 1989| , 



Tomisaka et al. 1989| , |Basu fc Mouschovias 199^) .f] To describe the transition between quasi- 



static evolution by ambipolar diffusion and dynamical evolution by gravitational collapse, 
[Li fc Shu (1996)1 introduced the idea of a pivotal state, with the scale-free, magnet ost at ic, 
density distribution approaching p oc for an isothermal equation of state (EOS) when 
the mass-to-fiux ratio has a spatially constant value, a condition that Shu fc Li (1997)| and 



Li fc Shu (1997)1 termed "isopedic" . Numerical simulations of the contraction of magnetized 



clouds justify the assumption of a nearly constant mass-to-flux ratio in the pivotal core.0 



^For example, as measured either by the net accelerations or by the square of the inward 
flow speed divided by the sound speed, the ambipolar-diffusion models in Figures 3 and 6 
of [Ciolek fc Mouschovias (1994)| spend less than 0.1% of the total computed evolutionary 
time in states where even a single grid point is more than 10% out of mechanical balance 
with self-gravity, magnetic forces, and thermal pressure (see also Figs. 3 and 7 of Basu & 
Mouschovias 1994 and Fig. 1 of Ciolek & Koenigl 1998). 

^For example, inside the starred point where [Ciolek fc Mouschovias (1994)| consider the 



core to begin, the mass-to-flux ratio varies in the last models of their Figures 3 and 6 by a 
factor of only 3 or 2 over a range where the density varies by a factor ~ 10^. Outside the 
starred point, the mass-to-flux value exhibits greater variation, but this occurs only because 
Ciolek fc Mouschovias (1994) impose starting values for the mass-to-flux in the envelope 



that are ~ 2 x 10 ^ times the critical value (see also Figs. 4a and 8b in Basu fc Mouschovias 



1994[ ). Such small ratios for the bulk of the mass of a molecular cloud are probably ruled 



out by the Zeeman OH measurements summarized by Crutcher (1998). 
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The small dense cores of molecular clouds that give rise to low-mass star formation are 
effectively isothermal ([Myers fc Benson 1983| ; ^hu, Adams, fc Lizano 1987| ). The situation 
may be different for larger regions that yield high-mass or clustered star formation. It has 
often been suggested that the EOS relating the gas pressure P to the mass density p of 
interstellar clouds can be represented by a polytropic relation P oc p^+^Z" with negative 
index n. [Shu et al. (1972}| pointed out the utility of this idealization within the context 
of the classic two-phase model of the diffuse interstellar medium [ [Fikerner (1967)| ; [b'ield,| 
Goldsmith, fc Habing (1969]| ; [Spitzer fc Scott (1969)|| , while [Viala fc Horedt (1974)| published 



extensive tables analyzing the stability of non-magnetized, self-gravitating spheres of such 
gases. pVlaloney (1988)| examined the linewidth-size and density-size relations of molecular 
clouds, first found by [Larson (1981)| and subsequently studied by many authors [e.g., [Leung 
et al. (19821 [Torrelles et al. (1983)[ , [Dame et al. (1986)[ , [Falgarone et al. (1992)[ , [Miesch 
fc Bally 1994|| . Maloney pronounced the results consistent with the properties of negative 
index polytropes. For a polytropic EOS, the sound speed Cg = {dP/dpY^'^ oc p^''^" increases 
with decreasing density if n < 0. The latter behaviour may be compared with the empirical 
linewidth-density relation for molecular clouds, Av oc p"'^, with q ~ 0.5 for low-mass cores 
( [Myers fc I'uller 1992| ) and q ~ 0.15 for high-mass cores ( [Caselli fc Myers 1995| ), implying 
that n lies between —1 and —3, or that a static F = 1 + 1/n lies between and 0.7. 

The case F = 1/2 is of particular relevance for the equilibrium properties of molecular 
clouds. [Walen (1944)[ found that the pressure of Alfven waves propagating in a stratified 
medium, Pwave oc |5Bp, in the absence of damping obeys the simple polytropic relation 
-Pwavc OC p^/^, a consequence of conservation of the wave energy flux t>yi|5Bp. This result was 



later derived more rigorously by [Weinberg (1962) in the WKB approximation for MHD waves 
propagating in mildly inhomogeneous media, and, more recently by [Fatuzzo fc Adams (1993) 



and [McKee fc Zweibel (1995)[ in a specific astrophysical context. In numerical simulations of 
the same problem, [Gammie fc Ostriker (199^ found indication of a much shallower relation 
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(r ^ 0.1) for a self-gravitating medium supported by nonlinear Alfven waves. On the other 
hand, for the adiabatic contraction of a cloud supported by linear Alfven waves, |McKee ^ 



Zweibel (1995)| found a dynamic 7 larger than 1. [Vazquez et al. (1997)| confirmed a similar 



behaviour in numerical simulations of the gravitational collapse of clouds with an initial field 
of hydrodynamic rather than hydromagnetic turbulence. 

In the limit of F ^ (or n — —1), the EOS becomes "logatropic," P oc In p, a 
form first used by P^izano fc Shu 19^ to mimic the nonthermal support in molecular clouds 



associated with the observed supersonic linewidths. The sound speed associated with the 
nonthermal contribution, Cg = {dP/dpY^'^ oc p~^/^ becomes important at the low densities 
characteristic of molecular cloud envelopes (as contrasted with the cloud cores) since the 
thermal contribution is independent of density if the temperature T remains constant. This 
nonthermal contribution decreases with increasing density and will become subsonic at high 
densities as recently observed in the central regions of dense cores ( Parranco fc Goodman 



[19981) . [McLaughlin fc Pudritz (1996)[ and [McLaughlin fc Pudritz (1997)[ have modeled the 



equilibrium and collapse properties of unmagnetized, self-gravitating, spheres with a pure 
logatropic EOS and claim to find good agreement with observations. 

Adams, Lizano, & Shu (in 1987) independently obtained the similarity solution for the 
gravitational collapse of an unmagnetized singular logatropic sphere (SLS), but they chose 
not to publish their findings until they had learned how to magnetize the configuration in 
a nontrivial way (see the reference to this work in Fuller & Myers (1992), who considered 
the practical modifications to the protostellar mass-infall rate introduced by "nonthermal" 
contributions to the support against self-gravity). Magnetization constitutes an important 
program to carry out if we try to justify a nonthermal EOS as the result of a superposition 
of propagating MHD waves (see also Holliman & McKee 1993). In this paper, we extend the 
study of Li & Shu (1996) to include the isopedic magnetization of pivotally self-gravitating 



- 6 - 



clouds with a polytropic equation of state. As a by-product of this investigation, we obtain 
the unanticipated and ironic result that the only way to magnetize a singular logatropic 
configuration and maintain a scale- free equilibrium is to do it trivially, i.e., by threading the 
SLS with straight and uniform field lines (see §6). 

A basic consequence of treating the turbulence as a scalar pressure, coequal to the 
thermal pressure except for satisfying a different EOS, is that we do not change the basic 
topology of the magnetic field. This assumption may require reassessment if MHD turbulence 
enables fast magnetic reconnection (|Vishniac fc Lazarian 1999D and allows the magnetic 



fields of highly fiattened cloud cores (Mestel & Strittmatter 1967, [Barker fc IVlestel 19961) 



or pseudodisks ( Galli fc Shu 1993b| ) to disconnect from their background. Recent MHD 



simulations carried out in multiple spatial dimensions (e.g., [Stone, Ostriker, fc Gammie 



T998t jMacLow et al 1998t Pstriker, Gammie, fc Stone 1999[ ; [Padoan fc Nordlund 1999[) 



find turbulence in strongly magnetized media to decay almost as fast as in unmagnetized 
media. Such decay may be responsible for accelerating molecular cloud core formation above 
simple ambipolar diffusion rates ([Nakano 1998| , [Myers fc Lazarian 1998[ , [Shu et al. 1999| ). 



Although this result also cautions against treating turbulence on an equal footing as thermal 
pressure, we attempt a simplified first analysis that includes magnetization to assess the 
resulting configurational changes when we adopt an alternative EOS for the pivotal state. In 
particular, different power-law dependences of the radial density profile translate immediately 
to different time dependences in the mass-infall rate for the subsequent inside-out collapse 
( [Cheng 19781 , McLaughlin fc Pudritz 1997). 



The paper is organized as follows. In §2 we formulate the equations of the scale-free 
problem and show that each solution depends only on the polytropic exponent F and a 
nondimensional parameter Hq related to the cloud's morphology. In §3 we present the 
numerical results. In §4, §5, and §6 we discuss the limiting form of the solutions. Finally, 
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in §7 we give our conclusions and discuss the possible implications of our results for star 
formation and the structure of giant molecular clouds. 



2. SELF-SIMILAR MAGNETOSTATIC EQUILIBRIUM EQUATIONS 



To begin, we generalize the singular polytropic sphere in the same way that [Li fc Shu 
(1996)1 generalized the singular isothermal sphere (SIS). In the absence of an external bound- 



ary pressure, the only place the pressure P enters in the equations of magnetostatic equilib- 
rium is through a gradient. Consider then the polytropic relation 

dP 



dp 



Kp 



1) 



By integrating equation (P we recover for F = 1 the isothermal EOS, P = Kp, (where K is 
the square of the isothermal sound speed) and for F = the logatropic EOS, P = K In p. 

We adopt axial symmetry in spherical coordinates and consider a poloidal magnetic field 
given by 

1 / $ \ 
^ = ^^^ ^^^<^ ' (2) 

where $(r, 6) is the magnetic flux. Force balance along field lines requires 



V + 



1 



F - 1 



Kp 



-(i-r) 



(3) 



where V is the gravitational potential and h{^) is the Bernoulli "constant" along the field 
line $ = constant. Poisson's equation now reads 



l_d_ 



d(^ dr 



dr 



+ 



d 



^2 sin^ e de 



d^ de 



de 



whereas force balance across field lines reads 

1 (d'^<^ 192$ coted^^ 



IGvr^r^ sin^ 6* \ dr'^ de"^ de 



dh 



AnGp; 
(4) 

(5) 
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We look for scale-free solutions of the above equations by nondimensionalizing and 
separating variables: 



PC \ i/(2-r) 



$ = 4 
dh 



'^3-2r^^4-3r\ i/(2-r) 



QT/2 j 

'23r-2^(^2-2r\ i/(4-3r) 



(6a) 
(6b) 

""V vri-r$2-r ) ' ^^^^ 

where Hq is a dimensionless constant that measures the deviation from a force free magnetic 
field, and R{9) and (j){9) are dimensionless functions of the polar angle 9^ These assump- 
tions imply that the equilibria will have spatially constant mass-to- flux ratios (see below). 
Substitution of equation ( |6cD into equations (H) and (|) yields 

1 d 



sin 9 d9 



R 



(2 - r)^ 



and 



sin' 9 



(2 - r)^ 



-CrHoR(f) 



_(2-r)/(4-3r) 



(7) 



(8) 



where a prime denotes differentiation with respect to 9, and 



Ar = 2'^(3-2r)/(4-3r){2-r)^ 



Br = 2-(i-n(8-5r)/(4-3r)(2-r)^ 



^ 2-r(i-r)/(4-3r)(2-r)_ 



(9a) 

(9b) 
(9c) 



^These definitions are not applicable for F = 4/3 or F = 2. 
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In particular, for Hq = 0, eq. gives the dimensionless density for the non- magnetized 
singular polytropic sphere 

R 



4-3r 



i/(2-r) 

(10) 



_(2-r)2 

whereas eq. implies $ = for < F < 1, in order to satisfy the boundary conditions 
eq. (0). In this case, the mass-to-fiux ratio A,, is infinite. However, for F = 0, eq. (|]) admits 
also the analytic solution of $ oc sin^ 9 corresponding to a straight and uniform field, while 
the density function is R{0) = 1. Therefore, a spherical logatropic scale free cloud can be 
magnetized with a uniform magnetic field of any strength, and any value of the spherical 
mass to flux ratio is allowed.0 

For arbitrary values of F and Hq the ordinary differential equations (ODEs) and (§) 
are to be integrated subject to the two-point boundary conditions (BCs): 



lim sin 6 



0, 



0(0) = 0, 0'(7r/2) = O, i?'(7r/2) = 0. (11) 

The first BC implies that there is no contribution from the polar axis to the mass inside a 
radius r. The second BC comes from the definition of magnetic flux, i.e. no trapped flux at 
the polar axis. The last two BCs imply no kinks at the midplane. 

The equilibria are characterized by: 



a 



the spherical mass-to-flux ratio, Q 



A, ^ 2nG'/'J^^ = 2(i-r)/(^-r) ( R{e) sin Odd, (12) 
<I)(r,7r/2) V4-3Fy 0(7r/2) io ^ ' 



this case, = 2fi^ = [20(7r/2)2]-i. 
^The standard mass-to-flux ratio A = 27rG'^/^M($)/$ is not defined for the polytropic 
scale free magnetized equilibria because the integral Jq^'^ R{6)(f){6)^^ sin 6d6 diverges since it 
can be shown that R{6 = 0) 7^ for F < 1. 
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where M(r) is the mass enclosed within a radius r; 

(6) the factor D by which the average density is enhanced over the non- magnetized value 
because of the extra support provided by magnetic fields, 



' 4 


- 3r ■ 


-i/(2-r) 


i(2 


-r)2_ 


Jo 



D 

which is equal to 1 if i^o = (see eq. [Il0| ); 
(c) the sound speed, 



R{9) sine de, 



(13) 



(14) 



and (d) the Alfven speed 



(i-r)/(2-r) 



;^i/(2-r) 



+ 



4-3r 
2-r 



R{e) sin' e 



Both the sound speed and the Alfven speed scale as r° for F = 1, and r^^^ for V 
other values of F, the exponent of r lies between these two values. 



(15) 
0; for 



It is also of interest to define the ratio /i^ of the square of the sound speed and the 
square of the Alfven speed, each weighted by the density, which is a physical quantity that 
can be compared with observations: 



2 ^ /o^ cjpsinede ^ 2-r/(2-r). 



lo^"^ '^\P sill (^d(^ 



fV2 

Jo 



/ sin 9d9 



(16) 



If Cs represents only the thermal sound speed, then the observational summary given by 
Fuller & Myers (1992) would imply that /i^ ~ 1 in the quiet low-mass cores of GMCs, whereas 
/i^ ~ 10~^ in their envelopes. If we include in Cg, however, the turbulent contribution, then 
the turbulent speed is likely to be sub-Alfvenic or marginally Alfvenic, and /i^ ^ 1 everywhere 
is probably a better characterization of realistic clouds. 



3. RESULTS 



To obtain an equilibrium configuration for given values of F and i^O) equations (|^ and 
(P) are integrated numerically. The integration is started at = using the expansions: 
= ao^^ + R = bo + 62e^(^-^)/(^-=^^) + . . ., with e = sin^, and 62 = [(4 - 3r)/2(l - 
r)]ArHQafl'^ ^^^^^ ^'"^6g~^. The values of and 60 are varied until the two BCs at = 7r/2 
(eq. 0), are satisfied. For flattened equilibria (see below) it is more convenient to start from 
6 = tt/2, where the BCs 0'(7r/2) = and R'{7i/2) = are imposed, and integrate toward 
6 = 0. The values of 0(7r/2) and R{tt/2) are then varied until a solution is found that satisfies 
the two BCs at 6^ = 0. 

Figure 1 shows the resulting flux and density functions (f){6) and R{0) computed for 
Ho = 0.5 and values of F between 0.2 and 1. We reproduce the results of [Li fc Shu (1996) 
for F = 1, which is the only case that obtains perfect toroids (i.e., R[6 = 0] = 0); models 
with F < 1 have nonzero density at the polar axis. Figure 2 shows the corresponding density 
contours and magnetic field lines. In the limit F ^ 0, independent of Hq as long as it is 
nonzero, the pivotal configuration becomes thin disks with an ever increasing magnetic field 
strength. Table 1 shows the spherical mass to flux ratio Xr, the overdensity parameter D, 
and the ratio of the square of the sound and Alfven speeds /i^. This table shows that, for 
fixed Hq, /i^ decreases as F decreases because the magnetic field becomes stronger. For the 
same reason D increases. In contrast, A^, goes through a minimum as F decreases. Figures 1 
and 2 demonstrate that for F — (the logatropic limit), Hq is not a measure of the strength 
of the magnetic fields since diverges as R{0) 6{tt/2 — 6) (see §5 below). 

For fixed F, a sequence from small Ho to large Hq progresses through configurations of 
increasing support by magnetic fields, as demonstrated explicitly for the isothermal case by 
Li fc Shu (1996)| . This behavior is illustrated here for the F = 1/2 case in Figure 3, which 
shows the density contours and magnetic field lines corresponding to values of Hq from 0.05 
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to 1.5. Table 2 shows the corresponding values of A^, -D, and /i^. For small Hq, the equilibria 
have nearly spherically symmetric isodensity contours and weak quasiuniform magnetic fields 
that provide little support against gravity. With increasing Hq, the pivotal configurations 
fiatten. The case Hq = 1.5 is already quite disklike: the pole to equator density contrast 
is R{7t/2)/R{0) ^ 10^. For a thin disk, the analysis of ^hu fc Li (1997]| demonstrates that 



magnetic tension provides virtually the sole means of horizontal support against self-gravity, 
with gas and magnetic pressures being important only for the vertical structure. In the 
limit of a completely fiattened disk {Hq oo), A,. ^ 1 independent of the detailed nature 
of the gas EOS (see next section). Table 2 shows the spherical mass to fiux ratio A^, the 
overdensity parameter D, and the ratio of the square of the sound and Alfven speeds /i^. 
Again D increases monotonically and /i^ decreases monotonically as the magnetic support 
increases with Ho, while A^ goes through a minimum and tends to 1 for large Hq. 

Since the mass-to-fiux ratio A^ is a fundamental quantity that will not change unless 
magnetic field is lost by ambipolar diffusion, in Figure 4 we consider sequences where A^ 
is held fixed, but F is varied. This Figure shows the locus of the set of equilibria with 
Xr = 0.95, 1, and 2 in the {Hq, F) plane. Equilibria with A^ < 1 are highly fiattened when 
F ^ even for small but fixed values of Hq (see §6). In fact, to obtain incompletely flattened 
clouds when one takes the limit F ^ 0, one also needs simultaneously to consider the limit 
Hq 0. Unfortunately, because both the density and the strength of the magnetic field 
at the midplane diverge as the equilibria become highly flattened, we are unable to follow 
numerically the limit F ^ to verify if these sequences of constant A^ < 1 will hook to a 
finite value in the Hq axis, or will loop to ifo = 0, consistent with our demonstration in §4 
that fiattened disks do not exist in the logatropic limit .| 



^As the equilibria fiatten due to either small F or large Hq, it becomes necessary to 
determine the constants of the expansions of R{0) and (j){6) near the origin with prohibitively 
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We speculate that the results for < 1 have the following physical interpretation. 
According to the theorem of Shu & Li (1997), only if A itself rather than A^ is less than unity, 
the magnetic field is strong enough overall to prevent the gravitational collapse of a highly 
flattened cloud. However, for moderate Hq and F when A^ < 1, even the singular equilibria 
are probably magnetically subcritical, since there can be little practical difference between 
the spherical mass-to-flux ratio A^ and the "true" mass-to-flux ratio A for highly flattened 
configurations. The latter is formally infinite when F < 1 only because the mass column goes 
to zero a little slower than the field column when we perform an integration along the central 
field line (see footnote 3). In this interpretation, subcritical scale-free clouds with A^ < 1 
and intermediate values of F can become highly flattened because magnetic tension supports 
them laterally against their self-gravity while the soft EOS does not provide much resistance 
in the direction along the field lines. The squeezing of the cloud toward the midplane is 
compounded by the confining pressure of bent magnetic field lines that exert pinch forces in 
the vertical direction. Both the magnetic tension and the vertical pinch of magnetic pressure 
disappear when the field lines unbend, as they must to maintain the scale-free equilibria in 
the limit F ^ (see below). As a consequence, logatropic configurations become spherical for 
any value of A^. We leave as an interesting problem for future elucidation the determination 
whether there is still a threshold in A^ below which the SLS, embedded with straight and 
uniform field lines, will not collapse dynamically. 

4. THE THIN DISK LIMIT {Hq > 1) 

In the limit Hq 1, the cloud fiattens to a thin disk for any F < 1. Dominant 
balance arguments applied to the two ODEs of the problem reveal the following asymptotic 



increasing accuracy. 
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behaviour:[| 



(4-3r)/(2-r) 



- f{0)Hl^ 

To the lowest order in Hq the equation of force balance along field lines (eq. 0) becomes: 

1 d 



;i7) 



[11 



sin 9 dO 

4-3r 



sin^(Arr 



{2-r)/(4-3r) ji _ _g-(2-r)^/ 



-(i-r) 



4-3r 



5^j2(i-r)/(4-3r) 



_(2-r)2 V 2-r 

valid over the interval < < 7r/2, plus the the integral constraint 



(19) 



(2-rp;„ - ^i^)'B,/;/..-r./(«n,„«.o, (20) 

obtained by integrating eq. (0) from 6 = to n, and applying the first BC (eq. 0) on the 
polar axis. 



The constant Rq is proportional to the surface density of the polytropic disks, given by 

(21) 



S(r) = lim / or sin 9d9 , ^ 



/ u-^ \ i/{2-r) 



which, for F = 1 gives S —>■ Hqo^ /jiGr, as found by |Li fc Shu (1996) 



Eq. (P) expressing force balance across field lines reduces to 

1 [/" - cot ef + ^{^ + 1)/] = -c^r^'-^^'^'-'^^RoKo - ^/2), 



sin'^ 



(22) 



where the parameter £ is defined by 



+ = 



2(4-3r)(i-r) 
(2-r)2 



(23) 



''These expansions are not valid for F = 1. See Li & Shu (1996) for the correct asymptotic 
expansion in this case. 
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This is equivalent to the equation for force free magnetic fields 

f -cotef + £{£ + !)/ = 0, (24) 
valid over the interval < ^ < 7r/2, plus the condition 

2f (vr/2) = CrRof{7r/2)-('-^^/(^''^\ (25) 

obtained integrating eq. (|^) from 7r/2 — eto7r/2 + e, and taking the limit e — > 0. For integer i, 

(-) 

solutions of eq. (|^) regular at ^ = are Gegenbauer polynomials of order i and index |, 
(see e.g. Abramowitz & Stegun 1965). In general, it can be shown (Chandrasekhar 1955) 
that any axisymmetric force free field, separable in spherical coordinates, can be expressed 
in terms of fundamental solutions whose radial dependence is given by a combination of 
Bessel functions of fractional order, and the angular dependence by Gegenbauer polynomials 
of index |. In our case, the choice of T determines a particular exponent of the power-law 
for the radial part of the flux function, and hence the corresponding value of i (non-integer, 
except for F = and 1). 

Therefore, the magnetic field is force free everywhere except at the midplane where 
p 7^ and the condition of force balance across field lines has to be satisfied. In the thin disk 
limit discussed here, the boundary condition 0'(7r/2) = is clearly not fullfilled: the kink 
of (p at the midplane provides the magnetic support against self-gravity on the midplane. 
Currents must exist in the disk to support these kinks. 

With the definitions 

eq. (|19D transforms into 

z" + cot Oz' + £{£ + l)z = y" + cot Qy + i{i + \)y, (26) 
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which has the solution 

z{e) = y{e) + q{e), 

where q{9) is a solution of the homogeneous equation 

q" + coteq' + i{i+l)q = 0. (27) 

Therefore, 



4-3r 



8(6) 

and the integral constraint eq. (EO) becomes 



g_Ar^^/2(i-r)/(4-3r) 



-i/(i-r) 



(28) 



7r/2 C2 _ 

q{e)smede = ^^^-^^i^. (29) 



The problem is thus reduced to the solution of the two homogeneous equations eq. ( pif ) 
and eq. (^) for the functions f{6) and q{6) which are determined up to an arbitrary constant. 
However, the two integral constraints that would have determined these latter constants 
(eqs. ^ ^9|), contain the additional unknown parameter Rq. The system of equations is 
closed by the requirement that 

lim Xr = 1. 

Ho^oo 

Substituting eq. (|T^ and eq. (|TBp in eq. (|T2|), one obtains 

lim A. = 2(i-r)/(2-r) ( luL) = 1, 

^^o-oo V4-3ry 2/(7r/2) 



I.e., 

2-r 



Ro = 2'/^'-'H^^)f{^m, (30) 



which gives the remaining condition. 



Eq. (24) and (pTQ can be solved numerically by starting the integration at = with the 
series expansions: q{e) = qo[l-\e{e+l)e^+. . .], and f{e) = fo {O^ - |[£(£ + 1) + |]^^ + •.•}, 
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where go and /o are arbitrary constants.0 The constants go; /o and Rq are then determined 



by the constraints expressed by eqs. pO|), (|25|),and 



Figure 5 shows the functions f{6) and s{6) obtained for F = 1/2 and increasing values 
of Hq from 0.4 to 1.5 compared with the asymptotic expressions computed here. Aheady for 
Hq = 1.5, the actual f{6) and s{9) are very close to the corresponding asymptotic functions 
eq. {\i7\} and eq. ([T8|) . Table 3 shows the value of the angle a of the magnetic field with the 
plane of the disk, the flux function / evaluated at = 7r/2 (indicative of the magnetic field 
stength), and the surface density parameter Rq, as functions of F. The angle a ranges from 
45° for the isothermal case F = 1 to 90° in the logatropic case F = 0. Correspondingly, the 
magnetic flux in the disk and the surface density both diverge as F ^ for any large but 
finite value of Hq. 



5. THE QUASI-SPHERICAL LIMIT {Hq < 1) 



For the isothermal case F = 1, [Li fc Shu (1996)] have shown how the SIS is recovered 
for Hq 1 from a family of toroids with zero density on the polar axis. For F 7^ 1, in the 
limit Hq 1, the asymptotic expansions are given by: 



m 



4-3F 



1 i/(2-r) 



(2 - F)^ 



(4-3r)/(3-2r) 



+ . . . 



^(^)^(4-3r)/2(3-2r) 



+ . . . . 



To the lowest order in Hq, eqs. (|^ and (|]) become: 

1 d 



sin 9 dd 



sin 9 



^^^-(2-r)/(4-3r)^. 



(2-ry 

4-3F 



-V 



^The two original BCs on the function R{9) are of little use here: the one at = reduces 
to the condition lim6i^o(l — F)^^ sin^^g' = 0, trivially satisfied; the second BC, R'{7t/2) = 
cannot be applied because of the (5-f unction at n/2. 
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(31) 



and 



1 



4-3r 



1 i/(2-r) 



[g"-coteg' + iii+l)g] = -Cr 



9 



(2-r)/(4-3r) 



(32) 



sin^^ 



(2-r)2 



The BC for the functions p{6) and g{6) are the same as those for R{6) and (j){6) in eq. (p!l|). 

Figure 6 shows the convergence of the solutions of the full set of equations (0) and 
(I) obtained for F = 1/2 and decreasing values of Hq from 0.4 to 0.05, to the asymptotic 
solutions obtained by integrating the equations above. Notice that p{0) < and piji /2) > 
showing that the sequence of equilibria with F = 1/2 originates from the corresponding 
unmagnetized spherical state (eq. |1^) by reducing the density on the pole and enhancing it 
on the equator. The same behaviour is found for any value of F in the range < F < 1. For 
F = 1, the function p{6) diverges at 6' = 0, indicating that this expansion is not appropriate 
in the isothermal case, as in the case Hq ^ 1. For the same reason, the expansion also fails 
for F = 0, since both g{6) and p{6) diverge on the equatorial plane. 

These flattened configurations are supported by magnetic and gas pressure against self- 
gravity. The intensity of the magnetic field can become very high even though Hq is small, 
because the latter parameter measures not the field strength but the deviations from a force 
free field (see eq. |6cD . 



We consider in this section the logatropic limit F 0. As anticipated in § 2, for F = 



threaded by a straight and uniform magnetic field. This solution represents the only possible 
scale-free isopedic configuration of equilibrium for a magnetized cloud with a logatropic EOS. 
To show this, we use the results of § 4 and § 5 for Hq ^ 1 and Hq 1 to find the limit of 



6. THE LOGATROPIC LIMIT (F ^ 0). 



eq. (0) and (^ admit the analytical solution R = 1 and $ oc sin 



^ 9 corresponding to a SLS 
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the equilibrium configurations for F ^ and fixed (small or large) values of Hq. 



In the limit Hq ^ 1, T ^ 0, analytic solutions to equations ( p4D and (|27| ) exist. 
The magnetic field tends to become uniform and straight, f{9) /(7r/2) sin^6', but /(7r/2) 



diverges, as shown in Table 3, and therefore s{6 = 7r/2) also diverges (see eq. ^). Eq. (^Tj) 
shows in this limit that the surface density S is independent of r, therefore, no pressure 
gradients can be exerted in the horizontal direction. The value S = [K 127x0)^ ^"^RqHq 
diverges as F ^ for any value of Hq because limr^o-Ro = oo (see Table 3). The magnetic 
fiux threading the disk, = 2^^^"^ RqHqt'^ sin^ 9, becomes infinite in order to keep the mass 
to fiux ratio A,, equal to 1. Therefore, the limiting configuration approaches a uniform disk 
with infinite surface density, threaded by an infinitely strong uniform and straight magnetic 
field. 



If we now examine the case Hq <C 1, in the limit F — 0, it is easy to show from eq. 
that the magnetic field tends to become uniform, g{9) g{n/2) sin^^^, but limr^o 5'(7r/2) = 
oo. Consequently, the density function p{9) also diverges in 9 = 7r/2, and the configuration 
again approaches a thin disk threaded by an uniform, infinitely strong, magnetic field. 

We conclude that scale-free logatropic clouds cannot exist as magnetostatic disks except 
in some limiting configuration. In the absence of such limits, the equilibria are spherical and 
can be magnetized only by straight and uniform field lines; i.e., the magnetic field is force- 
free and therefore given by Hq = 0. The inside-out gravitational collapse of such a SLS 
would still proceed self-similar ly as in the solution of McLaughlin & Pudritz (1997), but the 
frozen-in magnetic fields would yield a dependence with polar angle that eventually produces 
a pseudodisk (Galh & Shu 1993a,b; Allen & Shu 1998a). 
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7. SUMMARY AND DISCUSSION 

We have solved the scale-free equations of magnetostatic equilibrium of isopedic self- 
gravitating polytropic clouds to find pivotal states that represent the initial state for the 
onset of dynamical collapse, as first proposed by [Li fc Shu (1996)] for isothermal clouds. 



Compared to unmagnetized equilibria, the magnetized configurations are flattened because 
of magnetic support across field lines. The degree of this support is best represented by the 
ratio of the square of the sound to Alfven speeds /i^, or the overdensity parameter Z), since 
they are always monotonic functions of Hq and F. 

Configurations with F = 1 become highly flattened as the parameter Hq increases. When 
F < 1 (softer EOS) the equilibria get flattened even faster at the same values of Hq, since 
along field lines there is less support from a soft EOS than for a stiff one. However, it seems 
that in the logatropic limit flattened disks do not exist: the singular scale- free equilibria 
can only be spherical uniformly magnetized clouds. Figure 7 shows a schematic picture of 
the (i^o; r) plane indicating the topology of the solutions for scale free magnetized isopedic 
singular self-gravitating clouds. 

In self-gravitating clouds, the joint compression of matter and field is often expressed 
as producing an expected relationship: B oc p'', with different theorists expressing different 
preferences for the value of n (e.g., Mestel 1965, Fiedler & Mouschovias 1993). No local (i.e., 
point by point) relationship of the form B oc p'^ holds for the scale-free equilibria studied in 
this paper. However, if we average the magnetic field strength and mass density over ever 
larger spherical volumes centered on r = 0, we do recover such a relationship: (B) oc (p)*^, 
where angular brackets denote the result of such a spatial average and k = F/2. 

We may think of the result {B) oc (p)^^"^ as arising physically from a combination of two 
tendencies, (a) Slow cloud contraction in the absence of magnetic fields and rotation tends 
to keep roughly one Jeans mass inside every radius r, which yields (p) oc {cD/Gr"^, or (p) oc 
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j,-2/(2-r) ^^2^ ^ {py~^- (b) Slow cloud contraction in the absence of gas pressure tends to 
keep roughly one magnetic critical mass inside every radius r, which yields {B)/r{p) oc A,. 
= constant, or (B) oc r{p) oc r~^^^'^~^^ oc {p)^^'^ if gas pressure (thermal or turbulent) plays 
a comparable role to magnetic fields in cloud support. Notice that our reasoning does not 
rely on arguments of cloud geometry, e.g., whether cloud cores flatten dramatically or not 
as they contract; nor does it depend sensitively on the precise reason for core contraction, 
e.g., because of ambipolar diffusion or turbulent decay. 

Crutcher (1998) claims that the observational data are consistent with k, = 0.47 ±0.05. 
If we take Crutcher's conclusion at face value, we would interpret the observations as re- 
ferring mostly to regions where the EOS is close to being isothermal F fa 1, which is the 
approximation adopted by many theoretical studies that ignore the role of cloud turbulence. 
The result is not unexpected for low-mass cloud cores, but we would not naively have ex- 
pected this relationship for high-mass cores and cloud envelopes, where the importance of 
turbulent motions is much greater. Unfortunately, the observational data refer to differ- 
ent clouds rather than to different (spatially averaged) regions of the same cloud, so there 
is some ambiguity how to make the proper connection to different theoretical predictions. 
There may also be other mechanisms at work, e.g., perhaps a tendency for observations to 
select for regions of nearly constant Alfven speed, (va) oc {B) / {p)'^/'^ ^ constant (Bertoldi & 
McKee 1992). Thus, we would warn the reader against drawing premature conclusions about 
the effective EOS for molecular clouds, or the related degree to which observations can at 
present distinguish whether molecular clouds are magnetically supercritical or subcritical. 

If molecular clouds are magnetically supercritical, with greater than 1 by order unity 
(say, Ar = 2), then an appreciable fraction (say, 1/2) of their support against self-gravity has 
to come from turbulent or thermal pressure (Elmegrcen 1978, McKcc ct al. 1993, Crutcher 
1998). Modeled as scale- free equilibria, such clouds with p? of order unity are not highly 
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flattened (see Tables 1, 2 and Figs. 2, 3). Suppose we try gravitationally to extract a subunit 
from an unflattened massive molecular cloud, where the cloud as a whole is only somewhat 
supercritical, A,. ~ 2. If the subunit's linear size is smaller than the vertical dimension of 
the cloud by more than a factor of 2, which will be the case if we consider subunits of stellar 
mass scales, then this subunit will not itself be magnetically supercritical. Magnetically 
subcritical pieces of clouds cannot contract indefinitely without flux loss, so star formation 
in unflattened clouds, if they are not highly supercritical, needs to invoke some degree of 
ambipolar diffusion in order to produce small dense cores that can gravitationally separate 
from their surroundings. 

On the other hand, if molecular clouds are magnetically critical or subcritical, with 
Aj. < 1, then almost all scale- free equilibria are highly flattened, with appreciably less 
than unity. On a small scale, any subunit of this cloud, even subunits with vertical dimension 
comparable to the cloud as a whole, would also be magnetically critical or subcritical. For 
such a subunit to contract indefinitely, we would again need to invoke ambipolar diffusion 
to make a cloud core magnetically supercritical. Thus, although the decay of turbulence can 
accelerate the formation of cloud cores, the ultimate formation of stars from such cores may 
still need to rely on some magnetic flux loss (but perhaps not more than a factor of ~ 2) to 
trigger the evolution of the cores toward gravomagneto catastrophe and a pivotal state with 
a formally infinite central concentration. 



On the large scale, if GMCs are modeled as flattened isopedic sheets, |Shu fc Li (1997) 



proved that magnetic pressure and tension are proportional to the gas pressure and force 
of self-gravity. Their theorems hold independently of the detailed forms of the EOS or the 
surface density distribution. If GMCs are truly highly flattened - with typical dimensions, 
say, of 50 pc x 50 pc x a few pc or even less - then many aspects of their magnetohydrody- 
namic stability and evolution become amenable to a simplified analysis through the judicious 
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application and extension of the theorems proved by |Shu fc Li (1997]| (e.g., see Allen & Shu 
1998b, |Shu et al. 1999|) . This exciting possibility deserves further exploration. 
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Table 1. Parameters of Equilibria with i^o = 0.5 



T Xr D /i^ 
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0.00181 
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Table 2. Parameters of Equilibria with T = 1/2 
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Table 3. Parameters for the Thin Disk Limit 



r 
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/('r/2) 
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oo 


00 
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Fig. 1. — The flux and density functions (f){9) and R{9) computed for Hq — 0.5 and F = 1, 
0.8, 0.7, 0.6, 0.5, 0.4, 0.3, and 0.2. The isothermal case is indicated by thicklines. 

Fig. 2. — Density contours and magnetic field hues for Hq — 0.5 and F = 1, 0.8, 0.7, 0.6, 
0.5, 0.4, 0.3, and 0.2. The isodensity levels correspond to R{0) = 2^^, the isoflux contours 
correspond to (j){9) — (0.2A;)^, where /c = 0, 1, 2, 

Fig. 3. — Density contours and magnetic field lines corresponding to F = 1/2 and Hq = 
0.05, 0.1, 0.2, 0.4, 0.6, 0.8, 1.0, and 1.5. The isodensity and isoflux contours are the same as 
in Fig. 2. 

Fig. 4. — The locus of the set of equilibria with = 0.95, 1, and 2 in the {Hq, F) plane. 

Fig. 5.— The functions f{e) and s{9) for F = 1/2 and Hq = 0.4, 0.6, 0.8, 1.0, 1.2, and 1.5, 
compared with the asymptotic expressions obtained for 1, shown by the thick dashed 

lines. For clarity, the functions are shown over the range < ^ < tt. 

Fig. 6.— The functions g{e) and p{e) for F = 1/2 and Hq = 0.4, 0.3, 0.2, 0.1, and 0.05, 
compared with the asymptotic expressions obtained for Hq ^ 1, shown by the thick dashed 
hues. 

Fig. 7. — Schematic picture of the {Hq, F) plane indicating the topology of the solutions for 
scale free magnetized isopedic singular self-gravitating clouds. 
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